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Abstract 

A search is performed for Higgs bosons decaying into invisible final states, produced in 
association with a Z° boson in e + e~ collisions at energies between 183 and 209 GeV. The 
search is based on data samples collected by the OPAL detector at LEP corresponding to 
an integrated luminosity of about 660 pb~L The analysis aims to select events containing 
the hadronic decay products of the Z° boson and large missing momentum, as expected 
from Higgs boson decay into a pair of stable weakly interacting neutral particles, such 
as the lightest neutralino in the Minimal Supersymmetric Standard Model. The same 
analysis is applied to a search for nearly invisible Higgs boson cascade decays into stable 
weakly interacting neutral particles. No excess over the expected background from Stan- 
dard Model processes is observed. Limits on the production of invisibly decaying Higgs 
bosons produced in association with a Z° boson are derived. Assuming a branching ratio 
BR(h° — > invisible) = 1, a lower limit of 108.2 GeV is placed on the Higgs boson mass 
at the 95% confidence level. Limits on the production of nearly invisibly decaying Higgs 
bosons are also obtained. 



(Submitted to Physics Letters B) 



The OPAL Collaboration 



G. Abbiendi 2 , C. Ainsley 5 , P.F. Akesson 7 , G.Alexander 21 , G. Anagnostou 1 , K.J.Anderson 8 , 
S.Asai 22 , D.Axen 26 , I.Bailey 25 , E. Barberio 7 ' p , T. Barillari 31 , R.J. Barlow 15 , R.J.Batley 5 , 

P.Bechtle 24 , T.Behnke 24 , K.W.Bell 19 , P.J.Bell 1 , G.Bella 21 , A. Bellerive 6 , G.Benelli 4 , 
S.Bethke 31 , O. Biebel 30 , O.Boeriu 9 , P. Bock 10 , M. Boutemeur 30 , S.Braibant 2 , R.M.Brown 19 , 

H. J. Burckhart 7 , S. Campana 4 , P. Capiluppi 2 , R.K. Carnegie 6 , A. A. Carter 12 , J.R. Carter 5 , 
C.Y. Chang 16 , D.G. Charlton 1 , C. Ciocca 2 , A. Csilling 28 , M. Cuffiani 2 , S. Dado 20 , A. De Roeck 7 , 

E.A.De Wolf 7 ' 3 , K.Desch 24 , B.Dienes 29 , J.Dubbert 30 , E. Duchovni 23 , G.Duckeck 30 , 
LP. Duerdoth 15 , E. Etzion 21 , F. Fabbri 2 , P.Ferrari 7 , F.Fiedler 30 , I. Fleck 9 , M.Ford 15 , A. Frey 7 , 
P. Gagnon 11 , J.W. Gary 4 , C. Geich-Gimbel 3 , G. Giacomelli 2 , P. Giacomelli 2 , M. Giunta 4 , 
J.Goldberg 20 , E. Gross 23 , J. Grunhaus 21 , M. Gruwe 7 , A.Gupta 8 , C.Hajdu 28 , M. Hamann 24 , 
G.G.Hanson 4 , A.Harel 20 , M. Hauschild 7 , C.M.Hawkes 1 , R.Hawkings 7 , G.Herten 9 , 
R.D.Heuer 24 , J.C.Hill 5 , D. Horvath 28,c , P. Igo-Kemenes 10 , K.Ishii 22 , H. Jeremie 17 , 
P. Jovanovic 1 , T.R. Junk 6 '*, J. Kanzaki 22 '", D.Karlen 25 , K.Kawagoe 22 , T. Kawamoto 22 , 
R.K. Keeler 25 , R.G. Kellogg 16 , B.W.Kennedy 19 , S.Kluth 31 , T. Kobayashi 22 , M. Kobel 3 '*, 
S. Komamiya 22 , T.Kramer 24 , A. Krasznahorkay Jr. 29 ' 6 , P. Krieger 6 '', J. von Krogh 10 , T.Kuhl 24 , 
M.Kupper 23 , G.D. Lafferty 15 , H. Landsman 20 , D.Lanske 13 , D. Lellouch 23 , J.Letts , 
L.Levinson 23 , J.Lillich 9 , S.L.Lloyd 12 , F.K. Loebinger 15 , J. Lu 26,b , A.Ludwig 3 '*, J.Ludwig 9 , 
W. Mader 3 '*, S. Marcellini 2 , A.J.Martin 12 , T. Mashimo 22 , P. Mattig" 1 , J.McKenna 26 , 
R.A.McPherson 25 , F. Meijers 7 , W. Menges 24 , F.S. Merritt 8 , H.Mes 6 ' a , N.Meyer 24 , 
A. Michelini 2 , S.Mihara 22 , G. Mikenberg 23 , D.J.Miller 14 , W. Mohr 9 , T.Mori 22 , A. Mutter 9 , 
K.Nagai 12 , I. Nakamura 22 '^, H.Nanjo 22 , H.A.Neal 32 , SW. O'Neale 1 '*, A. Oh 7 , M.J. Oreglia 8 , 
S.Orito 22 '*, C.Pahl 31 , G. Pasztor 4 ' 3 , J.R. Pater 15 , J.E. Pilcher 8 , J. Pinfold 27 , D.E. Plane 7 , 
O.Pooth 13 , M.Przybycieh 7 -™, A. Quadt 31 , K. Rabbertz 7,r , C. Rembser 7 , P.Renkel 23 , 
J.M.Roney 25 , A.M.Rossi 2 , Y.Rozen 20 , K. Runge 9 , K.Sachs 6 , T.Saeki 22 , E.K.G. Sarkisyan 7 '- 7 , 

A. D. Schaile 30 , O.Schaile 30 , P. Scharff-Hansen 7 , J. Schieck 31 , T. Schorner-Sadenius 7 ' 2 , 
M. Schroder 7 , M. Schumacher 3 , R. Seuster 13 '^, T.G. Shears 7 - 71 , B.C. Shen 4 , P. Sherwood 14 , 

A. Skuja 16 , A.M. Smith 7 , R. Sobie 25 , S. Soldner-Rembold 15 , F. Spano 8 '*, A. Stahl 13 , D. Strom 18 , 
R. Strohmer 30 , S.Tarem 20 , M. Tasevsky 7,d , R. Teuscher 8 , M.A.Thomson 5 , E. Torrence 18 , 
D.Toya 22 , I. Trigger 7 -"", Z. Trocsanyi 29 - 6 , E.Tsur 21 , M.F. Turner- Watson 1 , I.Ueda 22 , 

B. Ujvari 29 ' 6 , C.F. Vollmer 30 , P. Vannerem 9 , R. Vertesi 29 ' 6 , M. Verzocchi 16 , H.Voss 7 ' 9 , 
J. Vossebeld 7 ' fe , C.P.Ward 5 , D.R.Ward 5 , P.M. Watkins 1 , A.T.Watson 1 , N.K.Watson 1 , 

P.S.Wells 7 , T.Wengler 7 , N.Wermes 3 , GW. Wilson 15 -*, J.A.Wilson 1 , G.Wolf 23 , T.R.Wyatt 15 , 

S.Yamashita 22 , D. Zer-Zion 4 , L. Zivkovic 20 



1 School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 
2 Dipartimento di Fisica dell' Universita di Bologna and INFN, 1-40126 Bologna, Italy 
3 Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 
4 Department of Physics, University of California, Riverside CA 92521, USA 
5 Cavendish Laboratory, Cambridge CB3 OHE, UK 

6 Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ottawa, 
Ontario K1S 5B6, Canada 



1 



7 CERN, European Organisation for Nuclear Research, CH-1211 Geneva 23, Switzerland 
8 Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago IL 60637, 
USA 

9 Fakultat fiir Physik, Albert-Ludwigs-Universitat Freiburg, D- 79104 Freiburg, Germany 

10 Physikalisches Institut, Universitat Heidelberg, D-69120 Heidelberg, Germany 

11 Indiana University, Department of Physics, Bloomington IN 47405, USA 

12 Queen Mary and Westfield College, University of London, London El 4NS, UK 

13 Technische Hochschule Aachen, III Physikalisches Institut, Sommerfeldstrasse 26-28, D-52056 

Aachen, Germany 

14 University College London, London WC1E 6BT, UK 

15 School of Physics and Astronomy, Schuster Laboratory, The University of Manchester M13 
9PL, UK 

16 Department of Physics, University of Maryland, College Park, MD 20742, USA 
17 Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec H3C 3J7, 
Canada 

18 University of Oregon, Department of Physics, Eugene OR 97403, USA 
19 Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK 
20 Department of Physics, Technion-Israel Institute of Technology, Haifa 32000, Israel 
21 Department of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
22 International Centre for Elementary Particle Physics and Department of Physics, University 
of Tokyo, Tokyo 113-0033, and Kobe University, Kobe 657-8501, Japan 
23 Particle Physics Department, Weizmann Institute of Science, Rehovot 76100, Israel 
24 Universitat Hamburg/DESY, Institut fiir Experimentalphysik, Notkestrasse 85, D-22607 Ham- 
burg, Germany 

25 University of Victoria, Department of Physics, P O Box 3055, Victoria BC V8W 3P6, Canada 
26 University of British Columbia, Department of Physics, Vancouver BC V6T 1Z1, Canada 
27 University of Alberta, Department of Physics, Edmonton AB T6G 2J1, Canada 
28 Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P O Box 49, Hungary 
29 Institute of Nuclear Research, H-4001 Debrecen, P O Box 51, Hungary 
30 Ludwig-Maximilians-Universitat Munchen, Sektion Physik, Am Coulombwall 1, D-85748 
Garching, Germany 

31 Max-Planck-Institute fiir Physik, Fohringer Ring 6, D-80805 Munchen, Germany 
32 Yale University, Department of Physics, New Haven, CT 06520, USA 



a and at TRIUMF, Vancouver, Canada V6T 2A3 
b now at University of Alberta 

c and Institute of Nuclear Research, Debrecen, Hungary 

d now at Institute of Physics, Academy of Sciences of the Czech Republic 18221 Prague, Czech 
Republic 

6 and Department of Experimental Physics, University of Debrecen, Hungary 
f and MPI Munchen 

9 and Research Institute for Particle and Nuclear Physics, Budapest, Hungary 
h now at University of Liverpool, Dept of Physics, Liverpool L69 3BX, U.K. 
% now at Dept. Physics, University of Illinois at Urbana-Champaign, U.S.A. 



2 



J and The University of Manchester, M13 9PL, United Kingdom 

k now at University of Kansas, Dept of Physics and Astronomy, Lawrence, KS 66045, U.S.A. 

' now at University of Toronto, Dept of Physics, Toronto, Canada 

m current address Bergische Universitat, Wuppertal, Germany 

n now at University of Mining and Metallurgy, Cracow, Poland 

° now at University of California, San Diego, U.S.A. 

p now at The University of Melbourne, Victoria, Australia 

q now at IPHE Universite de Lausanne, CH-1015 Lausanne, Switzerland 

r now at IEKP Universitat Karlsruhe, Germany 

s now at University of Antwerpen, Physics Department, B-2610 Antwerpen, Belgium; supported 
by Interuniversity Attraction Poles Programme - Belgian Science Policy 

* now at Technische Universitat, Dresden, Germany 

u and High Energy Accelerator Research Organisation (KEK), Tsukuba, Ibaraki, Japan 

" now at University of Pennsylvania, Philadelphia, Pennsylvania, USA 

w now at TRIUMF, Vancouver, Canada 

x now at Columbia University 

y now at CERN 

2 now at DESY 

* Deceased 



1 Introduction 



The Higgs boson pQ is required by the Standard Model (SM) [2] but has not yet been ob- 
served [3J. At LEP II energies it should be produced mainly through the "Higgs-strahlung" 
process (e + e~ —>■ Z* — > H°Z°) if its mass is sufficiently low. In the SM, the Higgs boson dom- 
inantly decays into a pair of the heaviest kinematically accessible particles, which would be a 
b-quark pair at LEP II. In some models beyond the SM, however, the Higgs boson can decay 
predominantly into a pair of invisible particles if the process is kinematically allowed. 

The Minimal Supersymmetric Standard Model (MSSM) [I] is one of the models which allows 
for invisibly decaying Higgs bosons [5], through the h° — > XiXi process, where xi is the lightest 
neutralino, if the mass of xi is lighter than half of the Higgs mass and R-parity is conserved. If 
Xi is purely photino-like, the decay h° — ► XiXi is suppressed. In this case a decay h° — > XiX^ 
where x% * s the second lightest neutralino, becomes dominant if it is allowed kinematically. 
If the mass difference (AM) between X2 an d X? is small, the visible products of the decay 
X2 —> X?Z*/7 are soft and the event topology is similar to that produced by an invisible Higgs 
decay h° — ► XiXi- The h° — > X1X2 processes are therefore referred to as nearly invisible Higgs 
decays. 

In a non-linear supersymmetric model, the Higgs boson can decay into a neutrino plus a 
Goldstino [6], and the invisible decay can be dominant. In other models beyond the SM with 
a spontaneously broken global symmetry, the Higgs boson could decay into a pair of massless 
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Goldstone bosons, called Majorons [7], which couple strongly to the Higgs boson. In models 
with extra dimensions, the Higgs boson can decay into a neutrino pair [B] or can oscillate into 
invisible states if the Higgs boson mixes with a graviscalar [9] which is a scalar graviton and 
escapes in the extra dimension. Models which introduce hidden scalar sectors which couple to 
the Higgs sector can also cause invisible decays of the Higgs boson [10J. 

In this paper, a search for invisibly decaying Higgs bosons (h° — > is presented using 

the data collected at various centre-of-mass energies (y/s) between 183 and 209 GeV by the 
OPAL detector at LEP, corresponding to an integrated luminosity of 659.3 pb _1 . The topology 
of events containing invisibly decaying Higgs bosons produced through the process e + e~ — > h°Z° 
is characterised by the decay products of the associated Z° boson plus large missing momentum 
and a visible mass (M vis ) of the event consistent with mip. Here it is also assumed that the 
decay width of the Higgs boson is negligibly small. A search for invisibly decaying Higgs bosons 
with large decay width is presented in Ref. [UJ. The search presented here looks for a hadronic 
decay of the Z° boson in association with missing energy. To cover other Z° decay modes, the 
results from this search are combined with the results of the decay-mode independent h°Z° 
search [12] where the Z° decays into e + e~ or yU + /i~. The results of the invisibly decaying Higgs 
boson search at LEP I [13] are also included to enhance the sensitivities to lower Higgs boson 
masses. 

The same analysis is applied to search for the production of nearly invisibly decaying Higgs 
bosons: e + e~ — > Z°h° — > (qq)(x°X°') assuming a small mass difference AM = 2 and 4 GeV 
between x° an d X°' ■ The standard neutralino searches [H] are sensitive to cases with AM > 3 
GeV. Similar searches for an invisibly decaying Higgs boson have been carried out by the other 
LEP experiments |15j . 



2 The OPAL Detector, Data and Event Simulation 
2.1 The OPAL Detector 

The OPAL detector is described in detail in Ref. [TB]. The central tracking system consisted 
of a silicon micro-vertex detector, a vertex drift chamber, a jet chamber and z-chambers. In 
the range | cos^l < 0.73, 159 points could be measured in the jet chamber along each traclo 
At least 20 points on a track could be obtained over 96% of the full solid angle. The whole 
tracking system was located inside a 0.435 T axial magnetic field. A lead-glass electromagnetic 
calorimeter (ECAL) providing acceptance within | cos 6*| < 0.984, together with pre-samplers 
and time-of-flight scintillators, was located outside the magnet coil in the barrel region and at 
the front end of each endcap. The magnet return yoke was instrumented for hadron calorimetry 

a While motivated by the lightest neutralino of the MSSM, throughout this paper we use x° as a generic 
symbol for a neutral weakly interacting massive particle resulting from an invisible Higgs boson decay. 

b A right-handed coordinate system is adopted, where the a;-axis points to the centre of the LEP ring, 
and positive z is along the electron beam direction. The angles 9 and <fi are the polar and azimuthal angles, 
respectively. 
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(HCAL), giving a polar angle coverage of | cos#| < 0.99, and was surrounded by external muon 
chambers. The forward detectors (FD) and silicon-tungsten calorimeters (SW) located on both 
sides of the interaction point measured the luminosity and complete the geometrical acceptance 
down to 24 mrad in polar angle. The small gap between the endcap ECAL and FD was filled 
by an additional electromagnetic calorimeter, called the gamma-catcher (GC), and a counter 
consisting of tile scintillators called the MIP plug. 

2.2 Data and event simulation 



The search is performed using the OPAL data collected at i/i between 183 and 209 GeV with 
an integrated luminosity of 659.3 pb -1 . The integrated luminosities at each y/s are listed in 
Tabled 



Nominal ^1 (GeV) 


183 


189 


192 


196 


200 


202 


204 


205 


206 


208 


<V^> (GeV) 


182.7 


188.6 


191.6 


195.5 


199.5 


201.7 


203.7 


205.0 


206.5 


208.0 


Lumi. (pb v ) 


56.1 


178.2 


29.0 


71.7 


74.9 


39.3 


6.3 


71.4 


124.6 


7.8 



Table 1: Average centre-of-mass energies ((y/s) ) and integrated luminosities collected at each 
nominal centre-of-mass energy after detector status cuts. The uncertainty on the luminosity 
measurement is 0.5%. 



The signal detection efficiencies and expected number of background events are estimated 
using a variety of Monte Carlo (MC) samples. Signal samples of invisibly decaying and nearly 
invisibly decaying Higgs boson processes are produced using the HZHA generator [T7] . The 
Higgs bosons are produced in association with a Z° boson, and then are forced to decay into 
a pair of invisible particles. Samples of h° — > x°X° a t each \fs are produced in one GeV steps 
in the Higgs boson mass range from 1 to 120 GeV with 2000 events per mass point. The 
h° — > X°X°' samples are generated with mass differences AM of 2 and 4 GeV at 5 or 10 GeV 
Higgs boson mass intervals between 30 and 120 GeV, with 500 or 1000 events per point. The 
detection efficiencies are determined at fixed values of the Higgs boson mass using the above 
samples and then interpolated to arbitrary masses with a spline fit. 

The most important background processes are e + e~ — > W + W~ — > tvqq and e + e~ — > Z°Z° — > 
uuqq. The first of these channels fakes a signal when the lepton is within a jet or escapes 
detection along the beam axis, and the second is an irreducible background for Higgs bosons 
with masses in the vicinity of the Z° boson mass. The radiative multihadron process e + e~ — > 
qq(j) also contributes due to the escape of photon into the beam pipe. 

The background processes are simulated primarily by the following event generators. For 
two-fermion (2f) final states, events are generated by PYTHIA [TS] and KK2f [12] (qq(7)), BH- 
WIDE [22] and TEEGG [23] (e+e"^)), and KORALZ [21] and KK2f (//+//" (7) and t + t-( 7 )), 
for four-fermion (4f) final states, by grc4f [20] (4f processes with final states of e + e~ ff) and 
KORALW [21] (4f processes except final states with e + e~ //), and for so-called two-photon pro- 
cesses where the initial-state electron and positron radiate photons which interact to produce 
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additional final state fermions, by PHOJET [23], PYTHIA and Vermaseren [2HJ (hadronic and 
leptonic two-photon processes; e + e~qq and e + e~£ + £~). The generated partons are hadronised 
using JETSET [18] with parameters described in Ref. [27J. The resulting particles are processed 
through a full simulation [28] of the OPAL detector. 



3 Selection criteria 

The search criteria are optimised at each y/s using the MC samples with 10 mass points just 
below the kinematic limit for the invisibly decaying Higgs, h° — > x°X°- F° r the h° — > x°X°' final 
state, the decay products of the Z° may be accompanied by a soft jet with small visible mass 
and energy, aligned in the direction of the missing momentum. Since the two event topologies 
are very similar, the selection criteria for the h° — > x°X° are also applied to the h° — > x°X°' final 
states. The analysis begins with a preselection to ensure data quality, followed by a combination 
of cut-based and likelihood-based analysis. 

Experimental variables are calculated using the four-momenta of charged particle tracks, 
and ECAL and HCAL clusters. The clusters associated with tracks are also used in the energy 
and momentum calculations, after subtracting the momenta of tracks from the energy observed 
in the calorimeters to reduce double counting of energy [29|. 

3.1 Preselection 

The following requirements are applied to reduce beam-related background as well as two- 
photon events: 

(PI) The event must not contain any charged particle track or ECAL cluster with reconstructed 
energy greater than 1.3 x -Ebeam, where E beiLm is the beam energy. 

(P2) ^lis° S ^ V-^vis < 0.2, where E vis is the total visible energy and ^M^° s ^l >0 - 9 is the visible 
energy in the region defined by | cos#| > 0.9. 

(P3) N^ od /N ch > 0.2, where N%°* and N ch are the number of good charged particle tracks 
defined as in Ref. [30] and total number of tracks, respectively. 

(P4) M v ; s > 3 GeV, where M v i S is the invariant mass of the event. 

(P5) pt > 1.8 GeV, where px is the magnitude of the vector sum of the transverse momenta 
of the reconstructed objects in the event with respect to the beam direction. 

(P6) Forward energy veto: events are rejected if there is more than 2/2/5 GeV deposited in 
either side of the forward detectors, SW/FD/GC respectively, or if there is any significant 
activity in the MIP plug. This forward energy veto is introduced to ensure that the data 
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sample consists of well-measured events. The efficiency loss due to vetoes on random 
detector occupancy has been studied with a sample consisting of random triggers, and was 
found to be between 2.2% and 4.1%, depending on y/s. The signal detection efficiencies 
and the numbers of expected background events are corrected for such losses. 

(P7) iVj ets > 1, where iVj e ts is the number of jets reconstructed with the Durham algorithm pH] 
with a jet resolution parameter y cut = 0.005. This reduces monojet-like background caused 
primarily by beam-gas and beam-wall interactions. 

(P8) | cos^missl < 0.95, where # m i SS is the polar angle of the missing momentum of the event. 
The qq(7) background is reduced by this requirement. 

(P9) M^ iss > GeV 2 , where M^ iss is missing mass squared and is calculated with the visible 
mass scaled to the Z°-mass, i.e. M^ iss = s — 2^ 7 Ii z ° E v[s + m| . This formula is applied 

-'"vis 

to avoid a negative M^ iss . 

The number of data events remaining after these cuts and those expected from SM back- 
ground processes are summarised in the first row of Table [21 

3.2 Main selection criteria 

The main selection consists of a cut-based analysis followed by a likelihood-based analysis using 
the same technique as described in Ref. [32] • After the preselection (P1-P9) the following cuts 
are applied in sequence: 

(Bl) iVX d > 4. 
(B2) j9 T > 6 GeV. 

(B3) max(| cos 6^1) < 0.95, where 0j et is the polar angle of the jet axis after the event is forced 
into two jets with the Durham algorithm. This requirement leaves events containing well 
measured jets. 

(B4) The number of isolated charged leptons identified as in Ref [32] is required to be zero to 
reduce the background contribution from semi-leptonic W + W~ and Z°Z° events. 

(B5) 120 GeV > M vis > 50 GeV. 

The distributions of px and max(| cos#j et |), just before applying the respective cuts, are shown 
in Figure [TJ The numbers of selected events, the expected background and the signal efficiencies, 
after each cut, are shown in Table [21 

After applying the above cuts, the selected sample is divided into two categories, namely 
events with two jets ("2-jet") and with more than two jets (">2-jet"), where the number of jets 
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Cut 


Data 


Background 


Efficiency (%) 






Total 


2f 


4f 


m h o = 105 GeV 


Pre 


101653 


86767 


48277 


8299 


67.6 


Bl 


40031 


34158 


14253 


7513 


67.6 


B2 


16895 


17037 


10391 


6503 


65.8 


B3 


16694 


16882 


10372 


6379 


65.3 


B4 


11476 


11654 


8855 


2695 


56.1 


B5 


1045 


1069 


523 


532 


55.5 


LH 2-jet 


194 


205.7 ± 1.8 ± 1.3 


46.1 


157.1 


23.7 


>2-jet 


278 


279.4 ± 1.9 ± 1.5 


44.1 


233.5 


23.4 



Table 2: Cut flow table at ^Js =183-209 GeV. Each row shows the number of events after each 
cut of the selection (described in the text) for the data and the expected background. The back- 
grounds from two-fermion and four-fermion processes are shown separately. The contributions 
from two-photon processes are not shown individually but included in the total background. 
The background estimates are normalised to luminosity at each energy and summed. The first 
quoted error on background estimates is statistics and the second systematic. The last column 
shows the luminosity-weighted average of selection efficiencies for the Z°h° — > (qq)(x°X°) & na l 
state with m h o = 105 GeV. The last two rows show the final numbers of selected events, ex- 
pected background and the efficiency after the likelihood analysis (LH) in each category. The 
efficiency in a category is the fraction of signal Monte Carlo events which pass the selection 
requirements. The background numbers and signal efficiencies include the occupancy correction 
determined at each ^/s due to the forward energy veto. 
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is defined by the Durham algorithm with y cut = 0.005. A likelihood analysis (LH) is built up 
for each category separately, with the same technique as described in Ref. [22] using the input 
variables: 

• COS $miss 

• the acoplanarity angle when the event is forced into two jets, aC op- 

• the invariant mass of the two jets with the smallest opening angle, M^^f . 

• d,23, which is defined as E^ is xt/23, where 7/23 is the jet resolution parameter at the transition 
point from two to three jets in the jet reconstruction. 

• min(A^ t ), which is the smallest charged multiplicity of any jet in the event. 

The distributions of input variables for the expected background are different between the two 
categories as shown in Figure El due to the different contribution from background sources. 
The resulting likelihood distribution for each category is shown in Figure [31 The remaining 
background in the signal-like region is dominated by semi-leptonic 4f events. 

The properties of 4f background events are similar to the signal, thus broadening the like- 
lihood peak for the signal. The final results are obtained by requiring the likelihood to be 
larger than 0.2. The numbers of observed and expected events are summarised in Table El and 
Table El The efficiency in a category is defined as the ratio of the number of selected events in 
that category to the total number of produced events; the sum of efficiencies in the two cate- 
gories provides the total efficiency at a given mass point. The efficiencies for the h° — ► x°X 0> 
processes are relatively lower than those for the h° — > x°X°i as shown in Table HI 





(GeV) 


183 


189 


192 


196 


200 


202 


204 


205 


206 


208 


2-jet 


Data 


17 


52 


7 


19 


20 


12 


1 


22 


43 


1 


Background 


19.0 


54.6 


9.0 


22.4 


23.4 


12.2 


2.1 


22.1 


38.5 


2.4 


Eff.(%) 


17.7 


20.0 


21.1 


22.5 


27.2 


26.7 


28.5 


27.6 


27.2 


29.2 


>2-jet 


Data 


15 


78 


18 


31 


41 


12 


2 


32 


47 


2 


Background 


30.6 


76.1 


13.1 


31.7 


30.8 


15.4 


2.6 


28.5 


47.6 


3.1 


Eff.(%) 


17.5 


20.9 


20.7 


24.3 


25.0 


25.1 


24.9 


23.6 


26.1 


23.6 



Table 3: Number of candidate events and expected background for each category at each \/s, 
together with signal efficiencies for m^o = 105 GeV. 



The systematic errors on signal efficiencies and the numbers of expected background events 
are estimated using the following procedures. The uncertainty corresponding to the modelling 
of each selection variable is determined by comparing the mean values of the distribution of 
that variable between data and SM background MC samples at y/s = mzo after applying the 
preselection. Efficiencies and numbers of expected background events are estimated again, 
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10 


14 


17 


19 


22 


23 


24 


23 
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x u x u z* 
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12 


15 


18 


20 


23 


24 


24 


23 


24 


x u x u z* 
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12 


16 


19 


21 


23 


24 


25 
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Table 4: Luminosity-weighted averages of efficiency for x°X°> X°X°7 an d X°X°Z* channels. No 
Monte Carlo samples are available in the x°X°7 aR d X°X°Z* channels for Higgs masses less than 
30 GeV. 



shifting each variable separately by its uncertainty. Relative changes to the original values of 
efficiencies and numbers of expected background events are taken as systematic errors for that 
variable. The systematic errors for the LH selection are estimated in a similar way. The total 
systematic errors due to the modelling of the selection variables including those entering the 
LH selection are calculated by summing the errors in quadrature for each category at each ^/s 
individually. The evaluated errors are summarised in Table [5l The statistical errors due to 
the finite size of the MC samples and the uncertainty on the luminosity measurement are also 
estimated. The total systematic error ranges from 3.5% to 17.4% for the signal, and from 1.9% 
to 4.4% for the background. 



Category 


2-jet 


> 2-jet 


Decay Mode 


Inv. 


Nearly Inv. 


Inv. 


Nearly Inv. 


AM (GeV) 




2 


4 




2 


4 


Selection variable 


Signal 


0.2-3.9% 


0.0-5.0% 


0.0-11.1% 


0.3-4.1% 


0.4-4.7% 


0.6-10.3% 


BKG 


0.7-2.2% 


1.0-1.9% 


MC statistics 


Signal 


3.3-7.3% 


5.2-10.9% 


6.5-15.4% 


3.6-8.8% 


5.0-9.5% 


4.8-11.0% 


BKG 


1.9-3.8% 


1.4-2.7% 


Luminosity 


0.5% 


Total 


Signal 


3.5-7.8% 


5.3-10.7% 


5.8-17.4% 


3.7-9.5% 


5.1-7.8% 


4.8-11.4% 


BKG 


2.2-4.4% 


1.9-3.4% 



Table 5: Ranges of estimated relative systematic errors (in %) for all \fs. Errors on the signal 
efficiency are estimated at each MC mass point at each ^/i and those for background at each 
y/s. The total systematic error at each y/s is calculated by summing the individual errors in 
quadrature. 
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4 Results 



Figure H] shows the missing mass distribution for the selected candidate events together with 
the expected background and an expected signal of m h o = 105 GeV for the two categories, 
for all y/s combined. No significant excess above the expected SM background is observed in 
either category. The main background comes from four-fermion processes in both categories. 
The broad peak around 70 GeV in the four-fermion histograms is due to the e + e~ — > W + W~ 
process and the peak around 90 GeV is due to the e + e~ — > Z°Z° process. Final efficiencies are 
summarised in Table El 

Limits are calculated using the likelihood ratio method described in Ref. [33]. The M m i SS 
information is used as a discriminator in the calculation. Systematic errors on the background 
and signal estimate are taken into account. 



4.1 Limits on the production of invisibly decaying Higgs bosons 

Figure [5] (a) shows 95% confidence level (CL) limits on the production rate of an invisibly 
decaying Higgs boson relative to the predicted SM Higgs production rate, defined as 

a(e+e~ -> Z h°) 
BR(h° -> xV) 7+! 7 1 \ = BR(h° -> xVR 
a{e e — > h ri SM j 

where <r(e + e _ — > Z°h°) and <r(e + e~ — > Z°Hg M ) are the production cross-sections of the invis- 
ibly decaying Higgs boson and the SM Higgs boson, respectively, and BR(h° — > x°X°) is the 
branching ratio for the Higgs boson decay into a pair of invisible particles. 

The observed and expected ratios shown in the figure are obtained from the results of this 
search, combined with LEP I data [13] . and with results from the e + e~ and channels of 

the decay- mode independent searches [12]. For LEP I results, the recoil mass information is 
used as a discriminating variable, incorporated using a Gaussian mass resolution function; for 
the channels from the decay-mode independent search, the distribution of the squared recoil 
mass is used as a discriminant. 

The full line in Figure [5] (a) represents the observed upper limit at 95% CL on the relative 
production rate as a function of the Higgs boson mass. A Higgs boson which couples to the Z° 
boson with SM strength and which decays exclusively into invisible final states is excluded up 
to a mass of 108.2 GeV at 95% CL assuming BR(h° -> x°X°) = 100%, while a limit of 108.6 
GeV is expected. The compatibility of the data with the expected background is quantified 
using the confidence (p- value) for background-only hypothesis, 1 — CLb (see Ref. [3]) which is 
plotted in Figure E] (b). 
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4.2 Limits on the production of nearly invisibly decaying Higgs 
bosons 

The results obtained from two nearly invisible decay modes — > X°Z* and x°' ~^ X°l) i n 
this analysis are combined at each AM, where the lower of the two efficiencies is taken as the 
combined efficiency The limit calculation uses only results from this analysis. In Figure E] (a) 
and (b), limits on the production rate for a nearly invisibly decaying Higgs boson with AM = 
2 and 4 GeV are shown for the data taken between 183 GeV and 209 GeV. The production 
rate of the nearly invisibly decaying Higgs boson relative to the SM Higgs production rate is 
defined as 

BR(h° -> x°X°')R* 

where BR(h° — > x°X°') is the branching ratio for the decay into nearly invisible particles. The 
dependence of 1 — CLb on the Higgs mass for nearly invisibly decaying Higgs bosons is shown in 
Figure [6] (c) and (d). A Higgs boson coupling to the Z° boson with SM strength and decaying 
into the nearly invisible final states is excluded up to a mass of 108.4 and 107.0 GeV at 95% CL 
for AM = 2 and 4 GeV, respectively, assuming BR(h° — > x°X°') = 100%. The corresponding 
expected limits are 108.2 and 107.3 GeV. 



5 Conclusion 

A search for invisibly decaying Higgs bosons has been performed using the data collected by 
the OPAL experiment at centre-of-mass energies between 183 and 209 GeV, corresponding 
to an integrated luminosity of 659.3 pb _1 . The search has not shown any excess over the 
expected background from SM processes. Limits on the production of invisibly decaying Higgs 
bosons were calculated combining the results with those from LEP I and those from e + e~ and 
channels of a decay-mode independent search at LEP II. Invisibly decaying Higgs bosons 
with masses below 108.2 GeV are excluded at 95% CL if they are produced with SM cross- 
sections, assuming BR(h° — > x°X°) = 100%. The search criteria were also applied to a search 
for nearly invisibly decaying Higgs bosons. Limits of 108.4 and 107.0 GeV were obtained for 
AM = m x o/ — irtyQ = 2 and 4 GeV, respectively. 
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Figure 1: The distributions of pr & n d | cos#j Ct | at ^fs=20Q GeV before applying the sequential 
cuts on the variables: data (dots with error bar) are shown together with the predicted con- 
tributions from the background processes; e + e~ — > (cross-hatched), two-photon processes 
(negative slope hatched), four-fermion processes (positive slope hatched), and qq(7) (open). 
The background distributions have been normalised to 124.6 pb _1 . The distribution of sim- 
ulated signals for the process Z°h° — > (qq)(x°X°) f° r m h° =105 GeV are also shown with 
dashed line. The signal is normalised using the SM Higgs production cross-section and 100% 
production rate for the process Z°h° — > (qq)(x°X°)- 
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Figure 2: The distribution of likelihood input variables: cos6 miss , aC o P , , c?23 and 

mm(N 3 c h), at y/s = 206 GeV for the 2-jet category (left) and the >2-jet category (right). 
The background sources are shaded as in Figure H The distributions of the signal for simu- 
lated invisibly decaying Higgs bosons with m^o = 105 GeV are shown as dashed lines. The 
signal histograms are normalised to the number of events of the expected background. The 
first and last bins in each histogram include underflow and overflow, respectively. 
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Figure 3: The distribution of likelihood output for \fs= 183 - 209 GeV for (a) the 2- jet 
category and (b) the >2-jet category. The background sources are shaded as in FigurelH The 
distributions of the signal for simulated invisibly decaying Higgs bosons with m^o = 105 GeV are 
also shown. The signal histograms are normalised using 30 times the production cross-section 
of the SM Higgs boson and 100% production rate for the process Z°h° — > (qq)(x°X°). 




Figure 4: The distribution of missing mass for each category for all LEP2 data combined: 
(a) the 2- jet category and (b) the >2-jet category. The background sources are shaded as in 
Figure [21 The distributions of the signal for simulated invisibly decaying Higgs bosons with 
m\fi = 105 GeV are shown on top of the background distribution. The signal histograms are 
normalised using the production cross-section of the SM Higgs boson and 100% production rate 
for the process Z°h° — > (qq)(x°X°)- 
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m h o [GeV] 




m h o [GeV] 

Figure 5: (a) Observed and expected limits on the relative production rate for e + e~ — > Z°h° — > 
Z°X°X° (invisible decay) to the SM Higgs production rate at 95% CL as a function of the test 
mass m h o, assuming BR(h° — > x°X°) — 100%. The solid curves show the observed limits and 
the dot-dashed curves the median expected limits. The dashed and dotted curves show la and 
2o" bands of expected limits, (b) The background confidence 1 — CLb as a function of m\p. 
The thick solid curve shows the observed 1 — CLb and the thin solid curve the expectation in 
the signal plus background hypothesis. The dot-dashed, dashed and dotted lines show median 
1 — CL b , and the la and 2a bands expected for the background only hypothesis, respectively. 
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m h o [GeV] m h <> [GeV] 

Figure 6: Limits on the relative production rate for e + e~ — > Z°h° — > Z°x°X°' (nearly invisible 
decay) at the 95% CL, normalised to the SM production rate fore + e~ — > Z°H°, (a) for AM = 2 
GeV and (b) for AM = 4 GeV, assuming BR(h° — > x°X°') = 100% as a function ofm^o. 
Figure (c) and (d) show the 1 — CLb for AM = 2 and 4 GeV, respectively. 
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